KFKI-1994-10/A 
June 1994 



POLYNOMIAL FORM FACTORS IN THE 0(3) NONLINEAR a-MODEL 



J. Balog and T. Hauer 



Research Institute for Particle and Nuclear Physics 
H-1525 Budapest II4, P. 0. B. 49, Hungary 



ABSTRACT 



We study the general structure of Smirnov's axioms on form factors of local 
operators in integrable models. We find various consistency conditions that 
the form factor functions have to satisfy. For the special case of the 0(3) 
(T-model we construct simple polynomial solutions for the operators of the 
spin-field, current, energy-momentum tensor and topological charge density. 



In integrable quantum field theories physical quantities can be determined 
exactly in the bootstrap approach. Making general assumptions about the 
properties of the S-matrix has led to the complete determination of the fac- 
torized scattering matrix of many well-known models The bootstrap 

method has also been applied to the determination of matrix elements of 
local operators in these theories. This form factor bootstrap was initiated 
in [0-0 and further developed in refs. |]^-|T^. In addition to the explicit 



solutions for the Sine-Gordon, SU{2) Gross-Neveu and 0(3) a- models, the 
authors of refs. IQ'llOl ^-Iso gave a complete set of axioms that the form 



factors of an integrable model have to satisfy. Using these axioms (which we 
will call the Smirnov axioms), the form factor functions in other integrable 



models could also be calculated |]TT|-|jT3| 



In this letter we study the general structure of Smirnov's axioms and recast 
them into a form that is very convenient for the recursive determination of 
the form factor functions in the 0(3) NLS model. We shall show that apart 
from a simple multiplicative factor, which is explicitly given, the form factor 
functions in the 0(3) model are polynomials in the rapidity variables. These 
polynomials are determined by a simple recursive formula. Although we do 
not claim originality on the solution of the form factor functions we think 
that our approach will turn out to be useful for the practical study of the 
structure of operators in the 0(3) NLS model. (Note that while the solution 
for the form factor functions of the 0(3) NLS model is completely explicit 
in refs. [3-[|10|) the complicated definitions make them almost impossible to 



use in practical calculations.) 

For simplicity, we consider integrable models describing charge- self conjugate 
bosons of equal, nonzero mass m and no bound states, as in the 0(3) NLS 
model. (In case of bound states existing the analytical properties of the 
form factors are more complicated and model-dependent, and the existence 
of antiparticles makes the formulae hard to read because of upper and lower 
indices and charge conjugation matrices. We think, nevertheless, that it is 
possible to generalize the following arguments also for these cases.) 

We parametrize the asymptotic states by the rapidities and the internal in- 
dices of the particles. These states are created from the physical vacuum by 
the Zamolodchikov-Fadeev operators, whose exchange relation is governed 
by the (two-particle) S-matrix: 

z+(^)z+(^') = SABxxio - e')z+{e')z+{e) . (i) 

Here the capital indices belong to the internal degrees of freedom. 
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First we recall the list of requirements that in the bootstrap approach the 
two-particle S-matrix is assumed to satisfy. First, it is meromorphic in the 
complex 6'-plane and analytic in the < 96' < tt physical strip (no bound 
states). Bose-symmetry and T- reversal invariance imply: 

SAB,CD{d) = SBA,Dc{d) , (2) 

Sab,cd{G) = Scd,ab{.G) ■ (3) 

The most important assumptions in the S-matrix bootstrap are unitarity, 
crossing symmetry and Yang-Baxter equation, respectively: 

SAB,CDid)ScD,A'B'i — d) = ^AA'^BB' , (4) 
SAB,CD{iT^ - d) = SAD,CB{d) , (5) 

SAB,Vw{(^)SvC,A'zi(^ + (^')SwZ,B'C'i(^') (6) 

= Sbc,vw{G')Saw,zc'{0 + 6')Szv,A'B'{6) ■ 



Finally, we require that all the eigenvalues of the S-matrix at zero momentum 
transfer are —1 (which is satisfied in most integrable models, including NLS): 

5'ab,cd(0) = -5ad5bc ■ (7) 

Let us now take a local field operator X of spin s. We restrict our attention 
to its matrix elements between the vacuum and an n-particle state. (Other 
matrix elements can be derived from these ones [ll^.) The definition of the 
ra-particle form factor is: 

{o\x{o)\e,Au . . . , e^A^) = /S..^„(^i, . . . , . (8) 



The form factors are originally defined for ordered sets of real rapidities cor- 
responding to the asymptotic states, but they can be analytically continued 
to the complex plane in all variables. The Smirnov axioms postulate the 
properties of this analytically extended / function: 

• • • 5^n) = e^Vi"^..A„(^l - • • • , - A) , (9) 

f^%j...,e,e',...) = SAB,Yx{0-0')f3Y...{---,o',e,...), (lo) 

fiCL...AAOu02...,e^) = /il^^^^(e2,...,^.,^i-2^vr), (11) 

^^^/iBul.U.iP + /5, ^1, • • • , ^n) = ^ {SAB^tuAOu . . . , ^.) 

- Sbui...u„,Vi...v„a{6i, ■ ■ ■ ,()n\f3)fv^l,vA0i, . . . , 6'„)) . (12) 
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The matrix S{9i, . . . ,9n\l3) entering (|T2]) is a product of two-particle S- 
matrices corresponding to the scattering of particles {j3B,9iUi, . . . ,9nUn) 
into the set {9iVi, . . . ,9nVn, PA). One should also postulate that the func- 
tions / are meromorphic in all variables and they are analytic in the physical 
strip except for poles explicitly given in (|T^). (In case of bound states giving 
extra (nonkinematical) poles a fifth equation applies connecting the 1)- 
particle form factor to the n-particle one.) 

We can simplify eq. (|^) if we write / as a product of a "scalar" form factor F 
and an overall factor that carries the Lorentz-transformation character and 
cancels the e*^ in eq. (H): 

{0„...,9^)= (± e^») ' Fi^. {9,,..., 9^). (13) 



Equations (^^, (0) and (0) are of the same form (with F instead of /), 
except for (0) in the two-particle case, which changes to: 

(^^1 -O2- iT^yF%^^{9^, 92) = finite , {9^ ^9^ + ^vr) . (14) 



The first three axioms together describe an invariance property of the form 
factor function F, while (0) governs the structure of the poles. To identify 
the symmetry group described by the first three axioms, we summarize the 
way it is realized. Let us consider a group a manifold M., and the set of 
functions mapping the manifold into a linear space: T = {V : A4 C}. Let 
(f be an action of the group ^ on A^: 

^{v{P,9i),92) = ^{p,9i92) , \/gi,g2eG, peM (15) 

and assume that for all peAi, gtQ there corresponds a nonsingular linear 
operator Mg{p) acting on C satisfying 

Mg,{p)Mg,{^M) = Mg,g,{p) . (l6) 



Now we can define an action of the group on the set of linear functions, 
$ : ^ X JF ^ as: 

Mg,V)]{p) = Mg{p)V{^{p,g)), (17) 
$((?!, $((72, V^)) = ^{gig2.V), ygi,g2eg, peM . (18) 
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Introducing a compact notation for the action of (p and this can be sum- 
marized as: 



ip{p,g) = p^ (19) 

H9,V) ^ gV, (20) 

M,,ip)M,,ip'^) = M,,,M, (21) 

,V{p) = M,{p)V{p^). (22) 



The invariance group of the form factor functions can now be identified by 
studying the action on the arguments. (The elements of the manifold Ai 
are n-component rapidity vectors.) It is easy to see that the transformations 
form a direct product group with a continous and a discrete component. 
The continous component is the (l+l)-dimensional Lorentz group L acting 
through shifting all the rapidities by a common value as in eq. (^. The dis- 
crete component contains all the permutations of the n variables generated 
by the transpositions (p^Ol). In addition, any of the n rapidities can be shifted 
by 2ki7i, where k is an integer. (One can construct these transformations 
by repeatedly using (|T0]) and (0).) Since the permutations act on the gen- 
erators of this latter group as well, the discrete component is a half-direct 
product. (Note that this discrete group can be generated by the two elements 
corresponding to (|l^) and (0).) Thus the symmetry group of the first three 
axioms is: 

g = L(g){SnA Z^") . (23) 

It is easy to find the matrices corresponding to the linear transformations 
in eqs. (p|)-(pn]) explicitly and verify that as a consequence of the S-matrix 
properties (||)-(|^) they satisfy the consistency relation ([l6|). This means that 
the first three axioms can be summarized as: 

gF = F. (24) 

Now we turn to the fourth Smirnov axiom which determines the singularity 
structure of the form factors. We write the form factor F as a product of 
two terms, an overall factor which carries all the singularities and a regular 
term with the tensorial structure of F: 
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Here the function \1/ is chosen to satisfy 



Res{^iz);z = i7i) = , (26) 



7r2 



^{i7T + 0) = -^(ivr-^). (27) 

The advantage of using the reduced form factor G^"-* is that while it is still 
invariant under the group Q, it is analytic in the physical strip and the fourth 
axiom directly gives its value at the point 9i — 62 = in. More precisely, if we 
define 

^^^^ " mWT^)' ^^^^ 

RAB,CD{e) = -SAB,CD{0)E{e) , (29) 

SabM^) = RAB,cD{e)/E{-e), (30) 



then 5" plays the role of the "reduced" S-matrix and G satisfies (|T0D with 
S instead of 5* and (|TT|) with an extra (—1)"^^ factor inserted at the rhs. 
These modified equations, together with (^, define transformation matrices 
also satisfying the relations (|TB|) and they can be compactly written as 

^G = G, (31) 

where the upper left index notation now indicates this new, modified group 
action. In the rest of the paper we will use this modified action only. 



(|I2D now reduces to 

^ABC/3...C/„(/^ + ^^5/^5^3, ••• ,6'n) = TABU3■■■Ur^,V3■■■Vu{|3\(^■i■, ■ ■ ■ -i^^n) 

X 4::.l(^3,...,^n) (32) 

except for the n = 2 case which is modified to 

(^1 -02- i^y'^-'G^^l^M^ O2) = finite, (^i ^ 9^ + ^7r) . (33) 

The matrix appearing in (|3^) can be explicitly given in terms of the matrix 
R and the function E: 

TABU3...Un,V3...V„{P\03, ■ ■ ■,dn) = 

^{E{P -9s)... E{P - 9^)6abSu,v, • • • Su„v„ (34) 

+ (-1)"~^-RbC/3,X4V3(/5 _ 6)3) . . . RXr^U„,AVnif^ " ^n)) ' 
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Note that eq. ( p2D gives the value of the n-particle reduced form factor in 
terms of the (n — 2)-particle one at the specified point (and not its residue 
as (p!2D). Now (after choosing a suitable function our task is to determine 
a series of analytic functions G*^"^ that solves eqs. (l3ll) - (0). 

If we know the (n — 2)-particle reduced form factor, eq. (p2D gives the value 
of the n-particle one at the special point. Assuming that the n-particle form 
factor also solves the equations, one can apply the elements of the group Q 
and compute the values of this function at other special points as well. We 
can compute the value of G'-"-' at all points where any two of its arguments 
differ by an odd multiple of in. This opens the possibility of determining 
the form factors completely. If the form factors belong to some special class 
of functions, the knowledge of their values at some specially chosen points 
could be sufficient for their determination. This is the case for the 0(3) NLS 
model form factors, since the reduced form factors in this model, as we shall 
see, are polynomials. 

However, a question of consistency arises here. It is easy to see that the 
determination of the value of G^"-* at a special point (6*1, ... , 6n) defined above 
is not neccessarily well-defined, since in general there are many different 
points (^2 + ^2; • • • ; ^n) from which one can get there by applying the 
elements of ^. If the system of the equations is consistent, the different ways 
determining G*^"^ at a given special point should give the same result. This 
condition gives consistency equations. 

Now we present some equations of this kind, which will play an important 
role in the study of the 0(3)-model form factors. To make the formulae easier 
to read we introduce a compact notation: dot over (under) a variable means 
a shift by in {—in), and we denote the kih 6k rapidity simply by its index k. 
We define the following functions built from the {n — 2)-particle form factor: 

^Sc/3...[/„(12 ...n) = Tabc/3...c/„,Vi...V3(2|3 ...n) 

xGt5j3...n), (35) 
H(k) ^ P....P3^{2)^ fc = 3,...,n. (36) 

Here the PjS are those group elements that represent the transposition of the 
(z — l)th and the ith. component of the rapidity vector. The H^'^^s are nothing 
but the values of G*^"^ at special points: 

G^''\k23...n) = H^''\l2...n). (37) 
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We will use the following consistency equations satisfied by the H^'^^s: 







. 9k . . . 9i . . .) , dk — 9i , 


(38) 




= H^'^ 


s k, k + 1 , 


(39) 




= H^k-i) 




(40) 




= H'^k+i) 




(41) 


<)...^„(12...n) 


= (-ir 


'^£..A„A.(13...ri2), 


(42) 


<)..^Jl/3...n) 


= (-1)"^ 




(43) 



Using the identification (^) and assuming that G^"-* satisfies the Smirnov 
axioms it is easy to understand the meaning of these consistency equations. 
However, what we need is to prove them directly by using the definitions 
(|35D, ( p6D and the fact that the {n — 2)-particle form factors G^""^) entering 
these definitions satisfy the Smirnov axioms. Indeed, one can show that 
if the two-particle S-matrix possesses all the properties (0)-(§) and all the 
(n — 2) and the (n — 4)-particle reduced form factors (from which the if's are 
generated) satisfy the four Smirnov's axioms, then (|55|)-(^5D are satisfied. 

Before turning to the special case of the 0(3) model we note that the Smirnov 
axioms alone cannot determine the form factors completely. Indeed, it is easy 
to see that from a set of solutions G^^^ we can generate new ones if we multiply 
the form factors by scalar functions which are ^-symmetric 

((^1, . . . , 9ny) = fi(")(^i, ...,9n), ygeg (44) 

and satisfy 

f](«+2)(^ + ZTT, P,9,,..., 9n) = • • • , ^n) • (45) 

An important example of such functions is provided by the invariant squared 
mass of the n-particle state: 

/i2(^i, ...,9n) = (^cosh(^,))' - (5]sinh(^,))' . (46) 

From now on we consider the 0(3) a-model, where our considerations lead 
to the complete determination of the form factor functions. The S-matrix of 
the model is given by 0: 

Sab;Cd{0) = Si{9)6abScd + S2{9)6acSbd + S^{9)6ad6bc , (47) 
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where 

^2(0) = TTTT^^^TTT^^ ^ (48) 



2iTie 


{9 + i'K){9- 


2in) ' 


6(6 - in 


) 


{e + iTT){e- 


2z7r) ' 


2in{in — i 


0) 



{e + i'K){e - 2in) 
and we choose the well-known two-particle form factor for 

E{e) = {e + iTT){e - 2tn) . (so) 

The fact that makes the form factors explicitly calculable in this case is that 
the function E is identical to the denominator of the S-matrix (|48|) . From 
this it follows that the i?-matrix is (— Ix) the numerator of S, thus it is a 
polynomial of 6. Since both E and R are polynomials, the matrix T entering 
(|32D is also a polynomial of the rapidities. One can easily check that its degree 
is (2n — 5) in /3 and that it is quadratic in the other rapidities. From this it 
follows that if G^""^-* is a polynomial, the function H^"^^ is also a polynomial 
and it is straightforward to prove that all the H^^h are polynomials. 

More precisely the following statement holds: if G*^"^^) satisfies the axioms 
and it is a polynomial in all its variables of maximum degree (n — 4), then 
the H^'^h are polynomials of maximum degree (2n — 5) in 6k and (n — 2) in 
the other arguments. 

Now it is a natural assumption that G*^"^ is also a polynomial and we can 
write the following Ansatz: 

n n /in 

G(")(12...n) = 5:if«(12...n)n7T:yY' (51) 

k=2 1=2 v'^V 

where 

(kl) = {Ok - Oi) . (52) 



(|51|) is the main result of this paper. Together with the definitions ( p5D and 
( PUD it provides a recursive polynomial solution for all the form factors of the 
model. In writing down (pTI) we have made use of the fact that a polynomial 
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of degree (n — 2) (which we assume is in the variable ^i) is determined 
by its values at {n — 1) different points. This is precisely the information 
contained in (p^). 

It can be proven that ( [^ ) indeed defines a polynomial solution of the Smirnov 
axioms. First, one can show that if the {n — 2)-particle form factor satisfies 
Smirnov's axioms and it is a polynomial of maximal degree (n — 4) in all 
variables, then the function G^") constructed by (^11) is a polynomial in all 
its variables of maximal degree [n — 2). Here one has to use (38) and the 
counting of the degrees is straightforward. One can then check that eq. (p!0|) 
is satisfied by (|5TD by verifying it at {n — 1) points (which is sufficient since 
both sides are polynomials of degree (n — 2)). Here one uses (^)-(^). (p!l|) 
is also satisfied at these points as a consequence of (^2]) and (^31) . Finally, 



the last Smirnov axiom is satisfied by ( pT| ) by construction. 

To summarize, (0) provides us with a whole family of form factor functions 
corresponding to the matrix elements of the local field operator between the 
vacuum and an increasing number of particles, provided a starting element 
of this family (corresponding to n particles) is known, it is a polynomial 
satisfying Smirnov's axioms and its degree is not higher than (n — 2). 

Studying the one- and two-particle form factors of the model, we have found 
that the four basic operators (spin- field, current, energy-mo menum tensor, 
topological charge) all have matrix elements of the type described above. 

We define the reduced form factors of these four operators, respectively, as: 

fl...AAl---n) = ^oG^ZXi^---n), (53) 

fZ.AAl---n) = {te^'^)^oGt:ti^---n), (54) 

fZ.AAl---n) = (^±f'^){±e^'^)wt!2^^ (55) 
/a....a„(1...^) = (/i^-l)v&oGS:^l(l...r.), (56) 



where 



i<j 



Note that the operators are defined by their Lorentz- and isospin transforma- 
tion character. In case of the Lorentz- vector and (l,l)-tensor the conservation 
is included in the definition. The current and E-M tensor have non-vanishing 
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form factors for an even number of particles only, whereas the spin and topo- 
logical charge operators have odd particle form factors only. Since there is 
no invariant isovector, the family of form factors of the topological charge 
begins with the three-particle case and that is why the /i^ — 1 factor has been 
introduced here to cancel the unwanted poles of the three-particle matrix 
element. 

The starting elements of the respective families are: 

G^r^'^^O) ^ 6aA, (58) 
GtA^''{e^,e,) ^ e^A.A,, (59) 

gT~m \9^,9,) ^ _ (l^ _ e^^ ^AU. , (60) 

G'aiA2A3(^1'^2,6'3) = ^AiA2A:i- (61) 

It is easy to see that these functions satisfy the axioms. The form factors of 
the current and topological charge are really polynomials of the right degree 
and the recursion works automatically. In the case of the spin and E-M 
tensor this condition does not hold, but one can show that the Ansatz (pl|) 
does give the next element of the family (i.e. the function we get satisfies 
the axioms as well) and the three-particle spin and four-particle E-M tensor 
form factors are polynomials of degree 1 and 2, respectively. This means that 
we can use our recursive construction also for these operators. 

We think that the recursive formula (^TJ) will prove to be useful in the problem 
of studying the correlation functions in the 0(3) NLS model. The problem of 
constructing the correlation functions by summing over an infinite number of 



possible intermediate states is discussed in refs. |11],[13| for the case of inte- 
grable models with no internal quantum numbers. Some preliminary results 
in the case of the more difficult problem of the 0(3) model are discussed in 
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